Tuning the metamagnetism of an antiferromagnetic metal 

J. B. Staunton, 1 M. dos Santos Dias, 1 ' 2 J. Peace, 1 Z. Gercsi, 3 and K. G. Sandeman 3 

department of Physics, University of Warwick, Coventry CV4 7AL, U.K. 
2 Peter Griinberg Institut and Institute for Advanced Simulation, 
Forschungszentrum Jiilich and JARA, D-52425 Jiilich, Germany 
3 Department of Physics, Blackett Laboratory, 
Imperial College London, SW7 2AZ, U.K. 
(Dated: March 8, 2013) 

Abstract 

We describe a 'disordered local moment' (DLM) first-principles electronic structure theory which 
demonstrates that tricritical metamagnetism can arise in an antiferromagnetic metal due to the 
dependence of local moment interactions on the magnetisation state. Itinerant electrons can there- 
fore play a defining role in metamagnetism in the absence of large magnetic anisotropy. Our 
model is used to accurately predict the temperature dependence of the metamagnetic critical fields 
in CoMnSi-based alloys, explaining the sensitivity of metamagnetism to Mn-Mn separations and 
compositional variations found previously. We thus provide a finite-temperature framework for 
modelling and predicting new metamagnets of interest in applications such as magnetic cooling. 



PACS numbers: 75.30.Kz, 75.10.Lp,75.50.Ee,75.30.Sg,75.80.+q 



The application of a magnetic field to an antiferromagnet can cause abrupt changes to its 
magnetic state [T] . While such metamagnetic transitions have been known and studied for a 
long time in materials such as FeCl2 [2J and MnF 2 [3], it is their association with technologies 
such as magnetic cooling [I] that has driven recent efforts to control metamagnetism in the 
room temperature range, and at accessible magnetic fields. 

In the particular case of magnetic cooling, a large, or even 'giant' magnetocaloric effect can 
arise at a first order metamagnetic phase transition such as is found in FeRh [5]. However 
hysteresis can be present at these transitions, resulting in significant inefficiency during 
magnetic cycling. Therefore, antiferromagnetic (AFM) metamagnets that sit on the border 
between first order and second order behavior are of great interest. They emulate the low 
hysteresis, high entropy change properties of two of the most advanced magnetic refrigerants, 
La-Fe-Si [6] and MnFe(P,As) [7], both of which are instead ferromagnetic (FM) metamagnets 
that have field-induced metamagnetic critical points. It is worth noting that tricritical points 
have also been examined in the context of other technologies, such as displays based on liquid 
crystals [8]. 

Metamagnetic transitions have previously been investigated principally in terms of either 
localised or itinerant electron spin effects. In antiferromagnetic insulators [SJ, magnetic 
field-driven phase transitions can be understood qualitatively using a localised (classical) 
spin Hamiltonian with simple pairwise isotropic exchange interactions, a source of magnetic 
anisotropy and a Zeeman term for the influence of an external magnetic field. For example, 
Nagamiya jH] showed that, if the localised spins of a helical antiferromagnet are pinned 
by anisotropy and crystal field effects to spiral around a particular direction, the effect of a 
magnetic field is to bring about a first order transition to a fan structure where the moments 
now oscillate about the field direction. At higher fields, the fan angle smoothly reduces to 
zero to establish a high magnetisation phase at a second order transition. However, if 
anisotropy effects are negligible, there is now no favored axis for the helical order and, as 
soon as a magnetic field is applied, the helix plane orients perpendicular to the field. As the 
field is increased the moments smoothly cant into a conical spiral towards the field's direction 
with a second order transition to a high magnetisation phase. Thus in this localised picture, 
first order metamagnetism relies on a source of anisotropy. 

Metamagnetism is accounted for differently in an itinerant electron system. Seminal work 
by Wohlfarth and Rhodes [10J, Moriya and Usami developed and extended by many 
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others, e.g. [121 EE], derived the coefficients of a Landau-Ginzburg expansion of the free 
energy of an antiferromagnet exposed to a uniform magnetic field in terms of the antifer- 
romagnetic order parameter, Am^, with wave-vector modulation q. The coefficients of the 
expansion are determined from a consideration of both Stoner particle-hole excitations pU] 
and spin fluctuations (TTJ [13] generated from the collective behavior of the interacting elec- 
trons [HI HH [15] . Field-induced tricriticality occurs when the quadratic and quartic terms 
in Arriq- both equal zero. More recently, from an analysis of a generic mean-field Hamiltonian 
describing a helical antiferromagnetic state in an applied field Vareogiannis [TB] has shown 
how an itinerant electron system can undergo a first-order 'spin-flip' transition. For a weak 
magnetic field the AFM polarization is parallel to the applied field and flips perpendicular 
only when the field exceeds a critical value [T7] . 

In this Letter we describe an ab-initio spin density functional theory (SDFT)-based 'local 
moment' theory for antiferromagnetic metals that blends these two complementary accounts 
of metamagnetism to suit materials where slowly varying 'local moments' can be identified 
from the complexity of the electronic behavior. A principal result is that, in the limit where 
magnetic anisotropy effects are small or neglected entirely, the local moments' antiferromag- 
netic order can still undergo a first order transition to a fan or ferromagnetic state owing 
to the feedback between the local moment and itinerant aspects of the electronic structure. 
We test our theory against a detailed experimental case study of CoMnSi-based tricritical 
metamagnets and show how it provides quantitative materials-specific guidance for tuning 
metamagnetic and associated technological properties. 

We start from a generalisation of SDFT [18] that describes the 'local moment' picture of 
metallic magnets at finite temperature [T9H22] . Its basic premise is a timescale separation 
between fast and slow electronic degrees of freedom so that 'local moments' are set up with 
slowly varying orientations, {e,}. The existence and behavior of these 'disordered local 
moments' (DLM) is established by the fast electronic motions and likewise their presence 
affects these motions. This mutual feedback can then lead to a dependence of the local 
moments' interactions on the type and extent of the long range magnetic order through the 
associated itinerant electronic structure. The DLM picture of the paramagnetic state can be 
mapped to an Ising picture with, on average, one half of the moments oriented one way and 
the rest antiparallel. However once the symmetry is broken so that there is a finite order 
parameter {rh{\ profile, e.g in a FM or AM state and/or when an external magnetic field is 
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applied, this simplicity is lost. Ensemble averages over the full range of non-collinear local 
moment orientational configurations, {ej}, are needed to determine the system's magnetic 
properties realistically [22J. We now develop this idea and extend the DLM theory to enable 
the effects of an external magnetic field to be modelled along with temperature and field 
dependence of the magnetic state. 

We consider a magnetic material in an external magnetic field B at a temperature T. The 
probability that the system's local moments are configured according to {e{\ is Pj({ej}) = 
exp[— (3Q({ii}, B)]/Z where the partition function Z = JT. J dij exp[— /3f2({ej}, B)] and 
P — {^bT)^ 1 . An ab-initio 'generalised' electronic grand potential Q({ei},B) is in principle 
available from SDFT [18] where the spin density is constrained to be orientated according 
to the local moment configuration {ii}. The thermodynamic free energy, which accounts 
for both the entropy associated with the orientational fluctuations of the local moments as 
well as that for the remaining electronic degrees of freedom, is given by F = —ksT\ogZ. 
Q({ei},B) thus plays the role of a local moment Hamiltonian but its genesis can give it 
a complicated form. Nonetheless by expanding about a suitable reference 'spin' Hamilto- 
nian Qo{£i} = J2i hi ' &i an d, using the Feynman Inequality [2_3J, we can find a mean field 
theoretical estimate of the free energy [18] 

F({m l },B,T) = (Q({e l },B)) { ^ } 
i J i 

— * 

where the probability of a moment pointing along on a site i is Pi(ii) = exp[— (3hi ■ 
ii}/ J exp[— (3h ■ ei]dii, so that the set of local order parameters, each of which can take 
values between and 1, {mj} = {J e,Pj(ej)(iej} = {(ej)}. Our DLM self-consistent mean- 
field theory of the statistical mechanics of the local moments can be seen as the natural 
counterpart of the DFT self-consistent description of the interacting electrons. The first term 
of Eq. [T]is the internal energy of the local moments, i.e. the average of the electronic grand 
potential over local moment configurations consistent with the order parameter profile {fhi}, 
the second is (-T) multiplied by the local moments' entropy and the last their interaction 
with a field, B. The sizes of the local moments, {yUj}, are determined self-consistently [TS] 
via the generalised SDFT. 
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The Weiss field at a site / is given by 

c _ d(n{{ei},§)) {mi} 



(2) 



<9m; 

To capture the itinerant electronic component of the problem coming from the overall spin- 
polarisation of the electronic structure, we approximate the hi as the first two terms of an 
expansion about a uniform distribution rh = 4 ^^rhi, i.e. 



d(Q) 



dm 



-E 



dm/dm- 
3 J 

i,j (m) • {rhj - rh) 



(rhj - fh) (3) 



"(2) 

where S\-(fh) is the direct correlation function tensor and describes effective interactions 
between the local moments that depend on the magnitude of long range magnetic order, rh. 

A solution of Eqs{l]to[3]at a fixed T and applied field B minimises the free energy, Eq. [lj 
Several solutions, {m|^},{mp^},- • • may be found and the one which produces the lowest free 
energy, Eq. [TJ describes the equilibrium state of the system, {fh^Equii.- Hence metamagnetic 
transitions can be tracked as functions of T and B - for a given T the solutions for increasing 
values of B can show a transition from, say, an antiferromagnetic to ferromagnetic state at a 
critical field, B c (e.g. see Fig. 1(b)). The material's spin-polarised electronic structure also 
depends on B and T and the state of magnetic order. 

We now follow through the general framework laid out in Eqsjl] to [3] for a putative 
helical antiferromagnetic metal. In order to elucidate the main aspects, we shall assume 
temporarily that magnetic anisotropy effects are small and can be neglected. In the absence 
of an applied magnetic field, the solution {rh{\ which produces the lowest free energy has 
the form rhi = Am^(cos(q ■ Ri)x + sin(<f • Ri)y). The helical axis, z, has no preferred 
direction in the crystal lattice. The order parameter Am 9 - increases from to 1 as T drops 
from T/v to K. On application of a magnetic field, B, defining, say, the x-axis of our 
coordinate frame, we can find numerical solutions of the mean field equations [T] to [3] of (i) 
distorted helical form: rhi = (m + Am,-) cos(<f • Rj)x + Am^sm(q ■ Ri)y, (ii) a conical helix: 
rhi = rnx + Amq-(cos(q ■ Ri)y + sin(q ■ R\)z), (iii) a fan state: rhi = m x + Am,j-cos(g • Ri)y 
and (iv) a ferromagnetic state with Am 5 - = 0. The free energies are Fdh(m, Am 5 -; B,T), 
F c h(m, Amf, B,T), Ff an (m,Amf,B,T) and FpM(m; B,T) respectively and we search for 
the lowest one at a given B and T. Consequently the relative difference between F dh and 



any of the others determines whether there is a first order metamagnetic transition or not 
and if so the temperature dependence of the critical field, B c . 

The local moment interactions (fa) of Eq. [i] set up the Weiss field and, along with 
the uniform component, h(fa), are the key elements of our theory. We determine these 
quantities ab-initio using relativistic, spin-polarised, multiple scattering (Korringa-Kohn- 
Rostoker, KKR) theory and the coherent potential approximation (CPA) to handle the 
electrons moving in fields set by the orient at ionally disordered local moments [22j [2U 125] . 
Here, for the first time, we account for the variation of S\-{m) with increasing extent of long- 
range magnetic order fa driven by spin-polarisation of the itinerant electrons which mediate 
the interaction between the local moments [26} [27]. Since spin-orbit coupling effects are also 
included [25] we can assess whether the neglect of magnetic anisotropic effects is reasonable. 
Details on how S\ - (fa) is calculated can be found in references [18], [25J EHJ [29] . 

Our previous experimental investigations have revealed a class of magnetic materials 
based on the orthorhombic CoMnSi metallic antiferromagnet to be an ideal testing ground 
for the theory [3"UH3"2"] . Those studies considered both the composition-dependent metamag- 
netism and pronounced magneto-elasticity in CoMnSi by extensive magnetic and structural 
measurements, including a characterisation of the anomalous temperature variation of struc- 
tural parameters in zero magnetic field. CoMnSi orders into an non-collinear, helical anti- 
ferromagnetic state in zero magnetic field at T/v ~380 K. In an applied field this transition 
becomes a metamagnetic one to a high magnetisation/ferromagnetic state at a temperature 
T t . As the applied field is increased T t decreases before going through a tricritical point 
at around 2 Tesla where an enhanced magnetocaloric effect is observed along with marked 
changes to the interatomic spacings. All three subsystems, Co, Mn and Si can be doped 
and substituted which also affects the magnetic properties [32]. This material is typical of 
many useful magnetic metals in that its magnetism possesses both localised and itinerant 
electron spin attributes [17] . Localised magnetic moments can be identified with the Mn 
sites whereas the magnetism associated with the Co sites is reflected in the long range spin 
polarisation of the electronic structure. Our measurements have also found that the mag- 
netism of CoMnSi is dominated by the behavior of the Mn moments. Their interactions, 
however, are delicately poised and depend on the spacing between them, the compositional 
environment in which they sit and the overall long range magnetisation. 

We previously modelled a commensurate approximation to the helical AFM state [31] 
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found in experiment and other calculations have examined an assumed FM state [33]. This 
Letter makes the first ab-initio exploration of the non-collinear helical magnetism of CoMnSi 
and its magnetic field and temperature variation. We start with self-consistent KKR-CPA 
calculations [3IH36] for the paramagnetic (m = 0) DLM state of CoMnSi. Local moments 
of magnitude \i ~ 3.0/i£ are established on the Mn sites which are very close to the mag- 
netisation per Mn site we find in calculations of CoMnSi in a FM state (m = 1) and we 
can assume that the sizes of the Mn local moments are rather insensitive to the orientations 
of moments surrounding them. We thus use the potentials within the 'frozen potential' 
approximation [2U EI] to study CoMnSi for probability distributions {Pj(ej)} with a set 
of values of the magnetic order parameter m,, each ranging between and 1. We use an 
angular momentum cut-off of I = 3, a muffin-tin approximation for the potentials and local 
spin density approximation for exchange and correlation. For our DLM calculations for a 
finite long range order parameter, m ^ 0, we consider a probability distribution on a fine 
grid of local moment orientations [221 [21] ■ No local moment forms on the cobalt sites in 
the DLM paramagnetic state whereas for finite m, a small magnetisation associated with 
each Co site is induced by the lining up of the Mn moments and the consequent overall spin 
polarisation of the electronic structure. For the FM state, m — 1, this is ~ 0.6/iB per Co 
site [SUES]. 

The onset of magnetic order in CoMnSi and the temperature at which this transi- 
tion takes place is obtained from an examination of S\-\m) for the paramagnetic state, 
m = 0. By calculating the lattice Fourier transform of this interaction, S^ 2 \q, 0) and find- 
ing the wave- vector, q ma x, where the largest eigenvalue, s^ 2 \q, 0), maximises, we find the 
magnetically ordered state that the system forms below our mean-field theory estimate of 
T/v = s ( - 2 \q max ,0)/3kB. If q ma x = a FM state is indicated whereas q max ^ points to an 
AF state. Fig. 1(a) shows s^ 2 '(q, 0) using structural data obtained by neutron diffraction [30J. 
These data show the Mn-Mn distances, di, to vary by more than 2% over the temperature 
range 100-400K. Our calculations show that CoMnSi should order into an incommensurate 
helical AFM state along the c-axis, set by the orthorhombic crystal structure, at T/v ~ 400K 
in good agreement with experiment (301 EI]- They confirm the expectation that spin-orbit 
coupling effects are small in CoMnSi, with only a very weak anisotropy favoring the Mn mo- 
ments to lie in the a-b plane. This agrees with the very small level of magnetic anisotropy 
seen in measurement |38j . 
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CoMnSi is evidently very close to a ferromagnetic instability as shown by the value of 
s ( 2 )(g = o, m = 0) being close to that of s^ 2 ^(g max ,m = 0), a convenient signature for a 
potentially useful metamagnet. Its marked magnetoelastic properties are confirmed by the 
tendency to order ferromagnetically taking over from the incommensurate ordering propen- 
sity as d\ is increased, which experiment also shows to happen with increase in temperature. 
When the bath of electrons in which the Mn local moments sit develops a long range spin 
polarisation as a magnetic field is applied and m increases, our calculations show that the 
magnetic interactions, S^ 2 \q,m), weaken significantly only strengthening again in the spin- 
wave limit of m ~ 1. The inset to Fig. 1(a) illustrates this. The magnetic anisotropy effects 
remain very small, ( e.g we find a uniaxial anisotropy K < O.OlmRy. per unit cell for the 
completely ferromagnetically ordered state with a small orbital moment of < 0.03/ib on the 
Mn sites). 
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FIG. 1: (a) The Mn local moment interactions, s^ 2 \q, rh = 0), for CoMnSi with structures measured 
in neutron diffraction experiments at four temperatures labelled by the Mn-Mn spacing, d\, in A. 
The four Mn moments within a unit cell tend to align with each other. The inset shows the 
variation of s^ 2 \q max ,rh) with FM order parameter m for the case where d\ = 3.07 A. (b) The 
variation of m versus applied field for several temperatures for a fixed lattice structure, in which 
d\ = 3.07 A. A tricritical point is indicated at 372K, 2 Tesla. 



The ferromagnetically-induced weakening of the magnetic interactions acts to promote 
the stability of the distorted helical magnetic state over other states when the magnetic 
field is applied [33] • There can therefore be a first order transition to a fan, conical helix or 
ferromagnetic state at critical field B c despite the absence of magnetic anisotropy. Fig. 1(b) 
illustrates this effect, and shows our solutions of the mean field theory equations [T] to [3] 
where we have used our calculated ab-initio uniform Weiss field h(m) [22] and 5^(g mm ,m) 
values. We point out that there are no adjustable parameters for these curves which show 
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a transition from a distorted helical AF state (small m, B) to a high magnetisation state 
above a critical field B c . The transition is first order up to a tricritical temperature of 372K 
and second order thereafter. 

Fig. 2 shows effect on the temperature dependence of the critical field, B c of either 
changing Mn nearest-neighbor separations, d±, (Fig. 2(a)) in CoMnSi or varying composition 
slightly away from CoMnSi (Fig. 2(b)). If the itinerant electron effect which gives S^ 2 \q, m) 
its m-dependence is neglected and only the effects of a weak on-site magnetocrystalline 
anisotropy are included, we find the values of B c to be 2 orders of magnitude smaller. B c 
increases sharply as <ii decreases. This is in line with measured pressure dependence of 
B c of CoMnSi [39]. B c is also very sensitive to compositional doping - the AFM tendency 
strengthens as electrons are removed so that B c increases, as shown by an application to 
CoMng 5 Cr 5 Si, whereas it weakens when electrons are added, as in the case of CogsNisMnSi. 
The CPA is used to describe the compositional doping on the Co and Mn sublattice as well 
as the local moment disorder on the Mn sites. 




Temperature (K) Temperature (K) 

FIG. 2: (a) The critical field B c for transition between a helical AF and ferromagnetic state 
versus T for CoMnSi for the four structures referred to in Fig.l. and (b) for CogsNisMnSi, and 
CoMn9 5 Cr 5 Si, with d x = 3.06 A. 

To test the theory presented here we have carried out further analysis of our experimental 
data to extract the temperature dependence of the critical magnetic fields of CoMnSi and 
its compositionally doped counterparts, Cog 5 Ni 5 MnSi and CoMng 8 Cr 2 Si. These results are 
shown in Fig. 3(a). The lower panel, Fig. 3(b) shows our theoretical estimates of the critical 
fields for these same three systems, where we have used the zero field experimental structural 
data. Overall the magnitudes and trends are given very well by the theory. The variation 
in temperature of the relative spacings between the Mn atoms is responsible for concave 
appearance of B C (T) in CoMnSi. The Ni-doped structure is shown to be extremely sensitive 
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FIG. 3: (a) The critical field B c versus T for CoMnSi, CogsNisMnSi and CoMnggC^Si measured in 
experiment, (b) B C {T) for the same three systems from theory, using the experimentally determined 
lattice structures [32J in zero field. The Ni- and Cr-doped systems are also shown using the CoMnSi 
structural data (dashed lines). 

to structure, as observed experimentally [32J. When B c > 2 Tesla, appears to be 
overestimated by the theory. This discrepancy could be partly mitigated by allowing for the 
magnetic-field induced increase in d\ which can inferred from our calculations of h{m) and 
S^(q max , fn) which show that ferromagnetic tendency strengthens with increase in d%. 

Detailed theoretical models are needed to use in concert with experimental studies if the 
behavior of highly sensitive magnetic materials is to be analysed and tuned. We have dis- 
cussed and demonstrated one such model here and showed how the temperature and field 
dependence of the metamagnetic transition of an antiferromagnetic metal can be affected by 
varying composition and atomic spacing. Useful magnetic metals typically have large mag- 
netic moments which are established locally. The Disordered Local Moment model works 
well for these [221 Ell EH] and we have set out a theory for metamagnetic transitions based 
on this picture. An important aspect for magnetic material design is to distinguish first 
and second order magnetic transitions and to get some control of the location of precious 
tricritical points. Our CoMnSi case study reveals how reducing the Mn-Mn separation, d\, 
increases B c . Small degrees of compositional doping which reduce the number of electrons 
have a similar effect. One notable feature in our modelling is the dependence of the interac- 
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tion between the local moments upon the extent of long range magnetic order. This comes 
from the change in the behavior of the itinerant electrons that mediate these interactions 
and produces a mechanism for first-order 'spin-flip' transitions even for cases where there is 
little magnetic anisotropy. This aspect has an important role in the analysis and design of 
adaptive, magnetic metals. 

We acknowledge financial support from the EPSRC (UK) (J.B.S, J.P), a FCT Portugal 
PhD grant SFRH/ BD/35738/2007 (M.d-S.D), the Royal Society (K.G.S) and EPSRC grant 
EP/G060940/1 (Z.G). K.G.S. thanks M. Avdeev and J. Bechhofer for useful discussions. 
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